Conductance fluctuations at the integer quantum Hall plateau transition 
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We study numerically conductance fluctuations near the integer quantum Hall effect plateau 
transition. The system is presumed to be in a mesoscopic regime, with phase coherence length 
comparable to the system size. We focus on a two-terminal conductance G for square samples, 
considering both periodic and open boundary conditions transverse to the current. At the plateau 
transition, G is broadly distributed, with a distribution function close to uniform on the interval 
between zero and one in units of e^ /h. Our results are consistent with a recent experiment by 
Cobden and Kogan on a mesoscopic quantum Hall effect sample. 
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I. INTRODUCTION 

One of the early surprises of the burgeoning field of 
mesoscopic physics some ten years ago was the observa- 
tion of Iftffie conductance fluctuations in small metallic 
samples.ElH Metals with phase coherence lengths exceed- 
ing their size were found to exhibit sample (or field) spe- 
cific fluctuations in their conductance. The magnitude 
of the fluctuations - of order jh - was essentially inde- 
pendent of the mean conductance, leading to the name 
"universal conductance fluctuations." Theoretical expla- 
nations are based on models of diffusiag electrons, in 
which localization effects can be ignoredfl Generally, this 
requires that the mean conductance is much larger than 
e^//i, a condition fulflUed in the experiments. 

One of the striking features of the plateau transitions 
in the quantum Hall effect, era is that the magnitude of the 
macroscopic longitudinal conductivity Oxx is both metal- 
lic - independent of temperature as T — > ~ and of 
order e^//i. Conventional localization effects are inoper- 
ative due to the strong applied magnetic field. For the 
transition from insulator to the first filled Landau level, 
the experimental valuegj for the conductivity tensor are 
consistent with 
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These macroscopic conductivities are self-averaged, since 
the sample sizes are much bigger than the phase coher- 
ence length. Several authors have sivf ft theoretical ar- 
guments in support of these values ,la'LrEl although it is 
unclear that the averaging process is appropriate to the 
experiment. 

Recently Cobden and Kogan have measured the con- 
ductance of a smaU-nuantum Hall effect sample, in the 
mesoscopic regimcEj They find large fluctuations in a 
two-terminal conductance near the plateau transitions, 
as they vary the carrier density with a gate voltage. 



Specifically, the conductance seems to be almost uni- 
formly distributed on the interval between zero and one 
in units of e^//i. In striking contrast to conventional 
metallic samples, the magnitude of the conductance fluc- 
tuations is comparable to the mean conductance, G ~ 
{l/2)e^/h. 

In this paper we compute the conductance fluctua- 
tions employing a simple network model of the inter-pger 
quantum Hall effect (IQHE) plateau transition.Einl3 We 
extract a mesoscopic two-terminal conductance and its 
sample (and fleld) specific fluctuations. Our results are 
entirely consistent with the Cobden and Kogan experi- 
ment. Right at the transition, the conductance distribu- 
tion function is essentially uniform on the interval from 
zero to one in units of /h. 

Our paper is organized as follows. In Sec. II we briefly 
review the network model, specifying the appropriate ge- 
ometry and boundary conditions. The results for the 
conductance and its distribution are presented in Sec. 
HI. Sec. IV is devoted to a brief discussion. 



II. THE NETWORK MODEL 



To model the IQHE plateau transition, m; employ 
Chalker and Coddington's network model.ll^ In this 
model the interactions between the electrons are ignored. 
In their original formulation, the impurity potential was 
assumed to be slowly varying on the scale of the mag- 
netic length. The semiclassical trajectories moving along 
equipotentials were modeled via ballistic chiral propaga- 
tion along the links of a network. Quantum tunneling 
at saddle points between nearby equipotentials was in- 
corporated via tunneling at node parameters, connecting 
two incoming and two outgoing links. For simplicity, the 
nodes and links were placed on a regular (square) lat- 
tice. Randomness was incorporated via phase factors for 
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propagation along the links, which were assumed to be 
independent and uniformly distributed between and 2tt. 

In this paper we focus on the behavior near the plateau 
transition. Being a continuous (second order) phase tran- 
sition, we expect that universal critical properties (in- 
cluding conductance fluctuations) should not depend on 
details of the model. Thus, for example, the results ob- 
tained should also apply to systems for which the poten- 
tial is not varying slowly on the scale of the magnetic 
length. Extensive numerical simulations that have ex- 
tracted the critical exponent v for the, diverging localiza- 
tion length support tiiis supposition.B For example, Lee, 
Wang, and Kivelsonta have shown that inclusion of ran- 
dom scattering at the nodes gives the same jialue for v as 
in Chalker's original random-phase model.113 Moreover, 
consistent estimates for v have also been obtained from 
other numerical approaches, such as Thoulfiss. number 
studies of lowest Landau level Hamiltonians.BO A more 
serious concern is the legitimacy of ignoring Coulomb in- 
teractions between the electrons. It is conceivable that 
interactions - particularly long-ranged Coulomb forces - 
might be a "relevant" perturbation at the noninteract- 
ing phase tMjasition (fixed point), leading to new critical 
propertiesJlj't£l However, experimental measurements of 
V seem to coincide with the non-interacting value. Be- 
ing practical, we ignore Coulomb interacticms, and adopt 
Chalker and Coddington's network model.E^ 

To be specific, we study a square lattice network of 
nodes and links, as depicted in Fig. 1. The network is 
connected to two leads - to the right and left. The dis- 
tance between the leads, measured in units of the network 
lattice spacing, is denoted by L^. Of interest is the two- 
terminal conductance between these two leads. The sam- 
ple has a width, denoted Ly, in the transverse direction. 
We consider two boundary conditions in the transverse 
direction: (i) periodic boundary conditions and (ii) open 
boundary conditions. The case of open boundary condi- 
tions corresponds closely to the experimental geometry of 
Cobden and KoganJl3 In this case, edge states dominate 
the transport in the IQHE plateau. 

Quantum tunneling at each node is represented by a 2 
by 2 matrix, 

/ Wout \ _ f cosh 6* sinh6' \ / Win \ 
\ Win J ~ \ sinhe* cosh 6* J \ Wout / ' 

where v and w represent complex amplitudes for incom- 
ing and outgoing electron waves to the right and left, 
respectively, of a given node (see Fig. 1). By construc- 
tion, this matrix conserves the current, \win\'^ + kinP = 
I'fi'outP + l^^outP- The node parameter 9 determines the 
degree of backscattering at the node. So, for example, an 
incident wave, say w-m, is back-scattered into Wout with 
probability tanh^ {9) . To make the model invariant under 
a 7r/2 spatial rotation, the node parameters take two val- 
ues, 9i and 92, in alternating columns, and are chosen to 
satisfy the condition sinh(6'i) sinh(^?2) — 1- Randomness 
is incorporated via random phase factors along links. 



From symmetry, the plateau transition occurs when 
tanh(6'i) = tanh(6'2), or 6* = 6*^ = ln(l + ^2). It is thus 
convenient to define a variable A, which measures the 
"distance" to the transition 

A = tanh(0i) - tanh(6l2). (2.1) 

This parameter lies in the range — 1 < A < 1, and van- 
ishes right at the plateau transition. 




FIG. 1. Schematic representation of the network model 
for a square sample with size = Ly = 4. The ar- 
rows indicate the direction of wave propagation along 
the links, and 6i and 62 specify scattering at the nodes. 
The two-terminal conductance G is measured between the 
right and left leads. 

For open boundary conditions in the transverse (y) di- 
rection, the nodes on the top and bottom edges are mod- 
ified to be 

( U-out ^ 1 W Win \ 

\win ) \Q I )\ Wout ) 

so that Wont = Vin and Win = Wout- In this case, the 
boundary breaks (lowest Landau level) particle-hole sym- 
metry, A ^ — A, just as the edges do in a real physical 
system. In the Hall plateau phase, A > 0, extended edge 
states confined to the top and bottom boundaries of the 
sample are expected. In the localized insulator, corre- 
sponding to A < 0, all states are localized, even near a 
boundary. In real systems, the presence of an edge state 
accounts naturally for the quantized Hall conductivity. 

With transverse periodic boundary conditions, all 
nodes (in each column) are identical. In this case, there 
are of course no edge states. The two-terminal conduc- 
tance is only sensitive to extended bulk states, present 
at the plateau transition. The finite system is particle- 
hole symmetric (A —A), as will be clear from the 
numerical results. 

We compute numerically the total transfer matrix T in 
the X direction, for a network of width Ly and length L^. 
This Ly by Ly matrix relates the incoming and outgoing 
amplitudes in one lead - the w's - to the amplitudes in 
the other lead - the w's, and can be written schemati- 
cally as, W — TV, where W and V are vectors with Ly 
elements. To extract the two-terminal conductance G, it 
is useful to write T in the form 

f Wont \ _ f A B\f Vin \ 
\Win J \C D )\Vont )' 
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where Via denotes an Ly/2 column vector of the incom- 
ing ampUtudes in the right lead, Vout the outgoing, and 
similarly for the other lead. This can be inverted to ob- 
tain the iS* matrix relating all the incoming to outgoing 
modes, as 

f Wont \ _ f BD-^ A - BD-^C \ ( VFin \ 

Kut ) \ -D-^c JyVi^j- 

The two-terminal conductance can be expressed in terms 
of the Ly/2 by Ly/2 transmission matrix t, which relates 
the incoming amplitudes of one lead to the outgoing am- 
plitudes in the other lead: Vout = tWin- Thus we have 
t = Z3^^.|-The two-terminal conductance follows readily 
from t as ,113 

G=— tr[tt+]. (2.2) 

Below we focus on the two-terminal conductance for 
square samples of size L = L^ — Ly—A, 8, 16, and 24. 
The distribution function is obtained by evaluating G for 
many different samples. We typically take 5 x 10* sam- 
ples. 



for the largest sizes, the same value as the experimen- 
tally measured macroscopic Uxx- But this is probably 
coincidental, since our averaging procedure is not appro- 
priate for macroscopic samples. With open boundary 
conditions, the mean conductance is not invariant under 
A ^ — A, due to the breaking of particle-hole symme- 
try by the boundaries. At the plateau transition, G is 
slightly larger than 1/2, roughly 0.65 for the larger sizes. 
The increase of G'(A = 0) when changing the boundary 
conditions from periodic to open can presumably be at- 
tributed to an additional contribution coming from edge 
currents. 
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III. RESULTS 

The mean conductance, denoted G, is obtained by av- 
eraging G over a large number of samples. In Fig. 2(a), 
G is plotted versus the control parameter A, for four 
different sample sizes, all with periodic boundary condi- 
tions in the transverse direction. (In this and later fig- 
ures, G is plotted in units /h.) As expected the mean 
conductance is largest at the plateau transition, A = 0, 
reflecting transport via bulk delocalized states, and falls 
off towards zero away from the transition. For the larger 
system sizes, the peak in G narrows. With open bound- 
ary conditions, the mean conductance rises from zero to 
one, as shown in Fig. 2(b). As the sample size increases, 
the "step" becomes less and less rounded. Together, the 
two sets of data for G, resemble experimental plots of 
the macroscopic conductivities, Uxx and Uxy, when plot- 
ted versus electron density. However, a direct comparison 
with macroscopic conductivities is delicate, since the en- 
semble averaging procedure that we are using may not 
be appropriate to the thermal averaging taking place in 
macroscopic experimental samples. 

From scaling arguments, one expects that the mean 
conductances, G{L,A), should be only a function of the 
single scaling parameter Jj L^/'^A for large enough L and 
small enough A. Here is the localization length critical 
exponent, which has a value v ~ 7/3. In Fig. 3, the 
the data for G in Fig. 2, are replotted versus the scal- 
ing variable L^^'^A, with v = 7/3. Although the smaller 
sizes show marked finite-size scaling corrections, the data 
collapse is satisfactory at the larger sizes. 

It is amusing that with periodic boundary conditions, 
G at the plateau transition A = 0, is very close to 1/2 



FIG. 2. The mean conductance G plotted vs A for 

square systems of four different sizes L = A{x), L — 8(A), 
L — 16(Q), and L — 24(0) with (a) periodic and (b) open 
boundary conditions. 
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L'/"A (1^ = 7/3) L'''"a (1^ = 7/3) 

FIG. 3. Scaling collapse of the mean conductance from 
Fig. 2, plotted vs the parameter AL^^" with u = 7/3, for 
(a) periodic and (b) open boundary conditions. 

To characterize the conductance fluctuations first con- 
sider the root-mean-square conductance, defined as 

SG=^CP-G\ (3.1) 

where the overbar denotes an ensemble average over dif- 
ferent samples. In Fig. 4, we plot 5G versus A for four 
different sample sizes, with periodic boundary conditions 
in (a) and open in (b). For both boundary conditions, 
there are large fluctuations near the plateau transition, 
with SG{A = 0) « 0.3. Away from the transition, the 
fluctuations drop off rapidly with increasing system size. 
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as expected. In both cases, the peaks sharpen with in- 
creasing L, as expected from finite size scaUng. The phys- 
ical origin of the shght double-peaked structure in Fig. 
4(a) is unclear. In contrast to the mean conductance 
itself, the root-mean-square conductance is relatively in- 
sensitive to finite-size effects at the transition. 




FIG. 4. Root mean square conductance 5G plotted vs 
A for a square system of four different sizes: L — 4(x), 
L = 8(A), L = 16(a), and L = 24(0) with (a) periodic 
and (b) open boundary conditions. 

More informative than the root-mean-square conduc- 
tance, is the full conductance distribution function, de- 
noted P{G). To obtain this, we simply make a histogram 
plot of the number of samples with conductance G, for a 
very large ensemble. In Fig. 5, we plot the conductance 
distribution function at the plateau transition (A = 0), 
for the largest system size L = 24, with (a) periodic and 
(b) open boundary conditions. In both cases, the conduc- 
tance is very broadly distributed, roughly uniform over 
the interval from zero to one /h. Above G = e^//i, the 
distribution functions drop off rapidly, although there is 
a slightly larger "tail" with open boundary conditions, 
presumably due to edge current contributions. 
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FIG. 5. Conductance distribution function right at 
the plateau transition, A = 0, for the system size L = 24 
with (a) periodic and (b) open boundary conditions. The 
distribution functions are normalized to unity. 



IV. DISCUSSION 



The conductance distribution functions obtained nu- 
merically compare favorably with those measured in the 
Cobden_and Kogan experiment on a mesoscopic Hall 
system.EJ In this experiment, the two-terminal conduc- 
tance of a small sample (0.6 x 0.6 ^vn?) was measured as 
a function of carrier density, by varying a gate potential 
Vg. Large fluctuations in the conductance were seen upon 
varying Vg through the plateau transitions, as shown in 
Fig. 6(a). In the Hall plateaus themselves, smaller fluc- 
tuations were observed. This behavior is consistent with 
our numerics for the root-mean-square fluctuations, 5G in 
Fig. 4, which are largest at the plateau transition. From 
the data with gate voltages near the plateau transition, 
Cobden and Kogan obtained a conductance distribution 
function, shown in Fig. 6(b), which is roughly uniform 
on the interval from zero to one. Under the "ergodic" as- 
sumption that varying the gate potential is equivalent to 
changing the impurity configuration (i.e., the sample), 
we can compare their distribution function with ours, 
which was obtained by taking an ensemble average at the 
plateau transition (see Fig. 5). The similarity is striking. 



(b) 
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FIG. 6. Cobden and Kogan's experimental data: (a) 
Two-terminal conductance plotted vs the gate voltage 
Vg. The thick solid line is the same data averaged over 
a Vg interval of 16 mV. (b) Conductance distribution of 
data points near the plateau transition, in the interval 



To more closely mimic the experimental procedure, we 
have computed the conductance for a given sample, as a 
function of A. This is shown in Fig. 7, for a square sys- 
tem of size L = 24 with both periodic and open boundary 
conditions. Notice the very large fluctuations in the tran- 
sition region. The behavior in Fig. 7b with open bound- 
ary conditions is very similar to the "raw" experimental 
data of conductance versus gate potential. 

In addition to extracting conductance fluctuations for 
square samples, we have studied systems with various dif- 
ferent aspect rations, Lx/ Ly. For aspect rations between 
roughly 1/3 and 3, the qualitative results are essentially 
unmodified. For very long sample, however, » Ly 
we start seeing effects of one-dimensional localization. 
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FIG. 7. Conductance plotted vs A of a single sam- 
ple for the system size L = 24 with (a) periodic and (b) 
open boundary conditions. The solid lines are the mean 
conductance, which was obtained by taking an ensemble 
average. 

In summary, it appears that the quantum Hall plateau 
transition provides an ideal arena for studying finite-size 
effects on random phase transitions. The agreement be- 
tween simple models of non-interacting electrons, and 
the experimental data is striking. Among the open is- 
sues is the role of Coulomb interactions, which have been 
ignored in the numerics. Will they change the critical 
behavior, and possibly modify the conductance fluctua- 
tions? Moreover, even without interactions, an analytic 
description of the transition is lacking, for either fluctua- 
tions or average properties. One can only hope that the 
relative experimental accessibility of the quantum Hall 
plateau transition, will spur further theoretical develop- 
ments. 
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